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Dual superconductivity is believed to be a promising mechanism for quark confinement and has been 
investigated on a lattice effectively by a particular gauge called the maximal Abelian (MA) gauge. 
We propose a new formulation of SU(3) Yang-Mills theory on a lattice based on a non-linear change 
of variables where the new field variables are expected to reduce to those of the Cho-Faddeev-Niemi- 
Shabanov decomposition in the continuum limit. By introducing a new variable, say color field, 
carrying the color direction with it, this formulation enables us to restore and maintain color symme- 
try that was lost in the conventional MA gauge due to the naive separation of the gauge potential into 
diagonal and off-diagonal components. An advantage of this formulation is that we can define gauge- 
invariant magnetic monopoles without relying on specific gauges to investigate quark confinement 
from the viewpoint of dual superconductivity. In this talk, we will present the relevant lattice formu- 
lation to realize the above advantages and preliminary results of numerical simulations to demonstrate 
the validity of this formulation. This SU(3) formulation is an extension of the SU(2) version already 
proposed by us in the previous conference. 
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1. Introduction 

Quark confinement is still an unsolved and challenging problem in theoretical particle physics. 
Dual-superconductivity is believed to be a promising mechanism for quark confinement in quantum 
chromodynamics (QCD).[|l|] Indeed, the relevant data supporting the validity of this picture have been 
accumulated by numerical simulations especially since 1990 and some of the theoretical predictions 
have been confirmed in the Maximal Abelian (MA) gauge; infrared Abelian dominance, magnetic 
monopole dominance and non-vanishing off-diagonal gluon mass, which are the most characteristic 
features for dual superconductivity. However, they are not yet confirmed in any other gauge than the 
MA gauge and the MA gauge breaks color symmetry. 

In this talk, we propose a new compact lattice formulation for the SU(3) Yang-Mills (YM) theory 
to establish the dual superconductivity picture for quark confinement in a gauge invariant way. This 
could be a lattice version of the non-linear change of variables (NLCV) for the YM gauge field in 
the continuum formulation (originally known as the Cho-Faddeev-Niemi-Shavanov (CFNS) decompo- 
sition). There the YM gauge field is decomposed into two parts, A^ = + in such a way 
that the "Abelian (diagonal)" part is dominantly responsible for the area decay law of the Wilson 
loop average, while the remaining "off-diagonal" part decuple in the low-energy (or long-distance) 
regime, thereby, leading to the infrared Abelian dominance. For performing non-perturbative studies, 
therefore, it is important to give a procedure of extracting such an "Abelian" part and the remaining 
part X^ from the original YM gauge field A^ also on a lattice. We construct the SU(3) lattice formu- 
lation by extending the SU(2) version proposed in our previous work ||, ^, ^] (For details, see [Q]). 
In the SU(2) case, we have succeeded to define two compact lattice variables Vx.^ and Xv,^ which play 
the similar role to the "Abelian" and "off-diagonal" parts in the continuum theory. These new variables 
enable us to define a gauge invariant magnetic monopole in the compact formulation which guarantees 
that the magnetic charge is integer- valued and obeys the Dirac quantization condition. Moreover, the 
infrared "Abelian" dominance and magnetic monopole dominance in the string tension were demon- 
strated by numerical simulations, together with the non-vanishing mass for the "off-diagonal" part. It 
is crucial to introduce a color vector field n(A') for maintain the color symmetry of the original YM 
theory. 



2. A new compact reformulation of SU(3) YM 

We construct an SU (3) lattice formulation by extending the SU (2) case [^, ^. Two color vector 
fields, Hx and m^, are introduced. They play a crucial role of maintaining the color symmetry of the 
original YM theory. A link variable Ux^n represents exponential of the line integral of a gauge potential 

along a link from xto x + ejX: 

Ux,ii = .^exp (^ig J ^ dx^X^{x)^ = exp(-/g£Ay_^), (2.1) 

where e denotes a lattice spacing and the path ordering operator. In explicitly estimating the naive 
continuum limit we adopt the midpoint (/ = x + ejX /2) definition for the link variable. Then we 
obtain an extended theory, M-YM {Ux,^, n^- ,mv), which has an enlarged gauge symmetry SU (3)(b x 
[SU (3)/?7 (See Figure |l| ). Under the gauge transformation 0^- = exp(/6jt), Q.x = exp{i(Ox) G 
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SU (3), the color fields n^, nix and Ux,^ transform as 

n.v ®nx = @xnx@l , ni;, ®m.v = ©xm^Ql , (2.2) 

f/v./l ^ ^x,fl — ^xUx,^^x+fi (2-3) 

New compact lattice variables Vx,^ and Xt.^i which correspond to the decomposed variables in contin- 
uum theory, and X^, should be represented by n^, nix and Ux.^ in a similar way to the SU (2) case. 
The lattice version of the decomposition condition is given by 

% [Vx,f,]nx = 0, [Vx,n]mx = 0, (2.4) 

Tr(n,X,,^ ) = 0, Tr (m,X,^ ) = 0, (2.5) 

where D^^[Vx,^]^x is a lattice version of the covariant derivative defined by 

D'^, [Vx,f, ] 0. : = ^ {Vx,^, 0.+^ - (l>xVx,^ ) . (2.6) 
Under the gauge transformation 

Vx,f, ^ = exVx,i,@l+^ , (2.7) 
Xx^^ '^Xx^^ = &xXxji&l, (2.8) 

the decomposition condition satisfies desired gauge transformations, D^^[Vx,^]<px —>■ Df^[®Vx^^]®<px = 
®x [Df,[Vx,^]^x] ©1+^, Tr(nA,^) = Tr(®n,®X,,^) and Tr(m,X,^) = Tr(®m,®X,,^). To define the 
equivalent theory in terms of the new variables to the original YM, the extended symmetry should be 
restricted to the same symmetry as the original YM, SU (3)0=9. For this purpose, we use a new MAG 
(nMAG) condition which is obtained by minimizing the functional: 

F„MAG[i^,0;t/.,M,n.,m,] =£{||D;^[XM]Vr||' + ||^)^Pt/x,M]®m4H (2.9) 

x.fl x,ll 



Since ( [2.4| ) and ( |2.5| ) are coupled matrix equations, it is difficult to obtain the general solution. 
Therefore, we consider a formula which reproduce the NLCV in the continuum theory. We adopt the 
midpoint definition for Vv,/j and the site definition for Xx^^ : 

Yx,p. = expi-igeYx'.^i) = ^txp {-ig j'^'^^ dx^'N ^{x)^ , (2.10) 
X,,^ =exp(-/^£X,-,^), (2.11) 

where Vx^^ could be the link variable represented by exponential of the line integral of like Ux^^i ■ In 
the naive continuum limit D^^^x,n\^x agrees with the continuum version up to 0{e^) : 



D'^[Vx,y]^x = d^^^-ig[Yx',^M + '^{dp^'^x'-ig^x',^M >V.',^} + 0(£2). (2.12) 



We take an ansatz; 

yx,ii = aUx,^+^inxUx,^+^mxUx,^+^?,Ux,^nx+^+^AUx.nmx+^ 

+ y\nxUx,nnx+^ + y2mxUx.nmx+^ + y^nxUx,iimx+^ + yAnxUx,nmx+ii, (2.13) 



Vx,^ = Px.li^n , Px,ii ■■= \jyx,iiVl^ , (2. 14) 

Xx,^:=Ux,^Vl^, (2.15) 
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Figure 1: (Left panel) The gauge symmetry of the new formulation. (Right panel) The correlation functions for 
the color vector field : < «^n^ >, where (AB) elemens (A < B) are plotted. 

where the variable V^,/! is represented by polynomials of [/^^^ , ni;,^ and n^^ and should satisfy the property 
of the gauge transformation (12), ( ^jf ) and (2^). Here ( 2.14 ) represents the polar decomposition 
for obtaining the unitary matrix Vx,^ from V^,^ , where P^^^ is a Hermitian matrix, = ^x.^i ■ By 



substituting ( |2.13D into ([2.4|), the coefficients of Vx,^ are determined as 

Vx,^i = aUx,n + ynxUx,^nx+^ + ymxUx,nmx+^ . (2. 16) 

Note that ( 2.13| ) is used instead of Vx^^ to determine the coefficients as a necessary condition, since the 



relation D^[yv.^]nv = is obtained if D^[Vv,^]ii.r = is satisfied. Then, by using ( 2.15 ) and (2^), the 
coefficients in ( [2.16 ) are determined as 7 = 6a up to 0(e) for Xy^^, and the overall factor a can be set 
to a = 1. In the continuum limit, Vy,^ is given by 



Vy,^ = Tr(Ay^^iiy)ny +Tr(Ay,^my)m.^-/ + \/g + l/g [5^my,my] , 



(2.17) 



which agrees with V^(x') in the continuum theory. can be defined in two ways; one is a definition 
on the midpoint, Xy ^ = Ay ^ — Vy^^, and the other X;i-_^ on a lattice site ixoraXx^ii in terms of Ay^and 

Xx,p, = exp(-/geX.v,^) = Ux,iiVl^ = exp f -ige (Ay,^ - Vy,^) - [Ay,^,Vy,^] +o(£^) j . 



3. Determining the new variables from lattice data 

In this section, we consider a procedure for obtaining Vx^p, and Xx^^ from numerical simulations. 
To calculate P^,^, in general, we need to diagonalize = Vx^jiVx^^i by solving an eigenvalue problem. 
Here we discuss the way to obtain Vx^^ using the nMAG condition. 

Suppose a gauge transformation 0^ G SU (3) diagonalizes the color vector fields Hx and vAx such 
that ©vHxQJ = and ©x^x^X = A*^. This is always possible since [nv,m,v] = is satisfied. Using 
this gauge transformation Ux,^ is transformed as ®Ux,p = @xUx,p®]c+p, and Vx^^ and Px^^ can also 
be diagonalized at the same time, since we obtain the relations [Px^jx,mx] = 0, [Px^fx,nx] = using 
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the relations Vx^^m^ = nijcVt^^ and V^^^n^ = nxVx^^i from (p!4|). Indeed, we can rewrite ( 2.16| ) to a 
diagonalized form 

®V,,^=0A,mC^.= [V-,M + 6A3X/.A3 + A^®[/,,^A«] =3[diag(®4^^ (3.1) 

where denotes an {AB) element of ^U^^^. Then, we have the diagonalized form of Px,^, and V^,^ 
is given by 

X. = % C = diag i^y^y^) ■ (3-2) 

Therefore, the decomposition of Ux,^ is reduced to the problem of finding out such as = @x- The 
solution &x can be obtained by the nMAG condition. Since FnMAG is invariant under the local gauge 
transformation of SU (3)(o=e, the solution of the nMAG condition selects a gauge orbit in the extended 
gauge symmetry. Therefore we can always choose such a gauge transformation that diagonalizes color 
fields, = A3 and ^m^ = Ag. Using this ©J, the nMAG condition is rewritten to the same expression 
as the conventional MAG condition(G = &x^)- 

FnMAC = Y,Tri^Ux,i,X'^Ux7^X^)+Y,Tr{^Ux,^X'^^Ux-^X^)+c.c. 

x,ll x.jl 

=i(r<'.i'+r«f.i'+r<'.i')- (3-3) 

x,ll 

Note that this nMAG condition does not fix the gauge of the original YM, ^Ux,^, but selects a gauge 
orbit along the local gauge symmetry SU{3)a)=e- The conventional MAG condition corresponds to 
a special gauge choice of ®nx = and ®mv = A^ on the gauge orbit. When we choose an overall 
gauge condition of the original YM theory, for example, the lattice Landau gauge, ^Ux,^ = Ux.^, the 
configurations of and mx are determined using the gauge transformation 0^: = Qx'^Ux^^&l^^, 

nx = @%@x, mx = @%&x, (3.4) 
/ G^^n Gj^22 Gj^33 \ _ 
Vx,, = e^diag T^yul^y ^ (3-5) 



4. Lattice data 



Numerical simulations are done using the standard Wilson action of SU(3) YM. The configurations 
are generated on a 16^ lattice at j3 = 5.7 using the Cabibbo-Marinari heatbath algorithm!^. After 
5000 thermalizing sweeps with the cold start, 120 configurations are stored every 100 sweeps. We 
choose the lattice Landau gauge (LLG) for the overall gauge fixing of the original YM theory. In 
gauge fixing procedure, we use the over-relaxation algorithm to update link variables by using the 
gauge transformation of SU(2) sub-groups in the SU(3) gauge transformation. In order to avoid the 
lattice Gribov copy problem in the both LLG and nMAG condition, we try to find out the configuration 
which absolutely minimizes the gauge fixing functional. In the process of minimizing the gauge fixing 
functional for Ux,^, we have prepared 16 replicas generated by random gauge transformations form 
Ux,^, and among them we have selected the configuration which have attained the least value of the 
functional. 
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histogram of M charge - 





Figure 2: (Left Panel) The histgram of monopole charges 2^'^ and 2^^^. The number of the monoples are 
plotted for 120 configurations. Each distribution with integer valued monopole charges are represented on the 
grids. This is a preliminary resut. (Right panel) The logarithm plot of correlation functions of the gauge potential 
and the new variables: Z)^, Dyy and Dxx- 



First, we check color symmetry of our new formulation, which is the global SU (3) symmetry 
to be preserved in LLC Under the global gauge transformation, the gauge fixing functional of LLG, 
PllgIs] = L.Y.^Tr(^?7v,^)> is invariant, while the color fields nx,vcix change their directions. Therefore, 
we measure the space-time average of the color vector fields = Tr(A'*ii;c) and = Tr{X^mx) and 
the correlation functions. The right panel of Figure |l| shows, for examples, the correlation functions 
for iiv. The lattice data show that the color symmetry is preserved: = 0, {m^) = and 

{nX) = S^^DNNil), (mim^) = S^^DMrnil), (^4mf) =0, (y = x + lfi, ^ = 4). 

Note that this preserving color symmetry is an advantage of our new formulation. 

Next, we define a gauge invariant magnetic monopole using the "Abelian" part V^./j in the similar 
way to the SU(2) case |^. Two kinds of the gauge invariant magnetic monopole currents (a = 1 , 2) are 
defined by 

k^^^:=le,,apd''e^:l, (4.1) 
^=^'-gTr((^/ + n. + ^m,) V,,„y,+„,py4„V^;^i) , (4.2) 

The gauge invariance of is clear by definition. Note that is the a-th element of the diagonal- 
ized expression of Vx^^Vx+^^yVi+y^^Vx^y , i.e., diag( exp{ig^e^@l,y), exp{ig^e'^ely), exp(ig^ e^@ly)) . 
The left panel of Figure § shows the histogram of the magnetic monopole charges, indicating that 
integer valued magnetic monopoles are obtained. 

Finally, we investigate the propagators of the new variables. The correlation functions (propaga- 
tors) of the new variables and the gauge potential of YM are defined by 

Dooix-y) := (O^W(D)^(j)> for = A^, ^,V:^, ^, X^,^, (4.4) 

where an operator 0^{x) is defined as a linear type, e.g., Ax'^^ = (^Ux,fi — t^vV) / i^^s) ■ The right 
panel of Figure ^ shows preliminary measurements of correlation functions of Daa j Dyv and Dxx ■ The 
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correlation Dw corresponding to the "Abelian" part dumps slowly and has almost the same dumping 
rate as Daa, while the Dxx corresponding to the "off-diagonal" part dumps quickly. This suggests that 
the "Abelian" part of the gluon propagator is dominated in the infrared region, the and mass generation 
of the "off-diagonal" gluon. 

5. Summary and discussion 

We have proposed a new compact lattice formulations of SU(3) YM theory as an extension of the 
SU(2) case [Q][|6|l[01- This formulation has enabled us to define a gauge invariant magnetic monopole 
in the compact formulation which guarantees that the magnetic charge is integer-valued and obeys the 
Dirac quantization condition. We have shown that the new variables can be obtained in any gauge of 
YM theory. It is crucial to introduce the color fields n(x) and m(x) to maintain the color symmetry 
of the original YM theory. We have performed the numerical simulations and measurements on a 
lattice. We have shown that color symmetry is preserved for the new variables, and the integer-valued 
gauge invariant magnetic monopoles are obtained. Though these results are preliminary, the lattice data 
suggest infrared "Abelian" dominance and mass generation of the "off-diagonal" gluon by investigating 
the propagators of the new variables. The mass generation of gluon can be investigated in the same 
way in the SU (2) case, and these are under investigation. 

Through these numerical simulations, we have shown that our new formulation enables the gauge 
independent investigation of the confinement mechanism to overcome the problems in the conventional 
studies based on a special gauge such as the MA gauge. Further studies such as "Abelian" dominance, 
monopole dominance and mass generation of gluons in LLG and also in the other gauges, are important 
to establish the dual superconductivity picture. 
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